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ABSTRACT 

Using the Owens Valley and Nobeyama Radio Observatory interferometers, 
we carried out an unbiased search for hot molecular cores and ultracompact 
(UC) Hii regions toward the high-mass star forming region G19.61-0.23. In 
addition, we performed 1.2 mm imaging with SIMBA, and retrieved 3.5 and 
2 cm images from the VLA archive data base. The newly obtained 3 mm image 
brings information on a cluster of high- mass (proto) stars located in the innermost 
and densest part of the parsec scale clump detected in the 1.2 mm continuum. 
We identify a total of 10 high-mass young stellar objects: one hot core (HC) and 
9 UC Hii regions, whose physical parameters are obtained from model fits to their 
continuum spectra. The ratio between the current and expected final radii of the 
UC Hii regions ranges from 0.3 to 0.9, which leaves the possibility that all 0-B 
stars formed simultaneously. Under the opposite assumption — namely that star 
formation occurred randomly — we estimate that HC lifetime is less than ~l/3 
of that of UC Hii regions on the basis of the source number ratio between them. 

Subject headings: Hii regions — ISM: evolution — ISM: individual (G19.61-0.23) 
— stars: early type — radio continuum: ISM 



1. Introduction 

It is widely accepted that high-mass stars (> 8Mq) form in clusters. However, identifica- 
tion of cluster members is complicated because zero-age main-sequence stars and high-mass 
protostars are deeply embedded in their molecular surroundings. It is hence impossible to 
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use a simple SED classification analogous to that applied to low-mass young stellar objects 
(YSOs). In this context, ultracompact (UC) 0.1 pc) Hii regions and hot cores (HCs) play 
a fundamental role (e.g., Kurtz et al. 2000). The former are very young Hii regions ionized 
by an 0-B star which is embedded in the natal cloud. The latter are dense, hot, compact 
molecular cores which host even younger (proto) stars, although their life time has not been 
established. A comparison between the estimated number of UC Hii regions in the Galaxy 
and the 0-B star formation rate yields a duration for the ultracompact phase of ~ 10^ years 
(Wood & Churchwell 1989 [hereafter WC89]) which is two orders of magnitude larger than 
expected from the classical expansion model for Hii regions (e.g., Spitzer 1978). Several 
models have been proposed to explain this discrepancy (e.g., Akeson & Carlstrom 1996, and 
references therein). Clearly, one must have better knowledge of properties of clusters to 
which HCs and UC Hn regions belong and of their natal molecular clumps. 

A well-studied example of high-mass (proto)star cluster is W49N (at a distance of 
11.4 kpc). Subsequent to the identification of 12 UC Hii regions powered by O stars (De 
Prcc ct al. 2000), Wilner et al. (2001) identified 6 HCs from an unbiased survey of line and 
1.4- mm continuum emission. Using the HC to UCHii region number ratio, they obtained 
an estimate of 25-100% for the lifetime ratio between the HC and UC Hii region phases. 
However, the W49N study has an important limitation that the region is very distant which 
sets hmits on the census of detectable YSOs. The high-mass star forming region G19.61- 
0.23 is an excellent target for a study of cluster formation because of its distance {d — 3.5 
kpc; Churchwell, Walmsley, & Cesaroni 1990), and potential richness in terms of YSOs as 
indicated by studies at cm (e.g., Garay, Reid, & Moran 1985 [hereafter GRM85]; WC89; 
Garay et al. 1998; Forster & Caswell 2000) and mid-infrared (MIR) (De Buizer et al. 2003) 
wavelengths. In this letter, we present a study of the cm contimmm emission with higher 
sensitivity and better angular resolution, complemented by mm wavelength observations. 



2. Observations and Archive Data Retrieval 

Aperture synthesis observations at 3 mm were carried out using Owens Valley Radio 
Observatory (OVRO) millimeter array and Nobcyama^ Millimeter Array (NMA) with a total 
of 6 array configurations in the period from 2002 December till 2004 April. More than 15 
molecular lines (R. Furuya et al. in preparation) were observed simultaneously with the 
continuum emission at 89.9 GHz. The flux density of the calibrator was bootstrapped from 
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planets with an uncertainty of 15%. Combining the continuum data from both arrays, we 
obtain a synthesized beam of 1'.'52 x l'/51, and an RMS noise level of 2.8 mJy beam~^. 

On 2003 July 29, we used the SIMBA bolometer array on the SEST 15-m telescope to 
obtain a parsec-scale map (3.0' x 4.0' in size) of 1.2 mm continuum emission. Observing rj 
Carinae and Uranus, we checked pointing accuracy (RMS ~ 2"), and made flux calibrations. 

We retrieved 3.5 and 2 cm data from the archive of the NRAO^ Very Large Array 
(VLA). The 3.5 cm data were taken from projects AW374A and AK470 (A-array), AK470 
(B), AK423 (CnB), and AK450 (D), and the 2 cm data from AW158 (B and C), AK423 
(CnB). Resulting beam sizes and noise levels are 0'.'33x0'.'23 and 0'/55x0'/50 with 0.31 and 
1.7 mJy beam~^ at 3.5 and 2 cm, respectively. 

3. Results and Analysis 

Figure 1 presents the SIMBA 1.2 mm continuum emission map where two clumps sep- 
arated by ~75" (i.e. 1.3 pc) are seen over a region of ~ 4.1 pcxl.7 pc. The brighter clump 
to the west hosts the G19.61-0.23 complex. A total flux density (Si,) of 19.4 Jy is obtained 
by integrating the emission inside the 3(7 contour level. Figure 2a shows the VLA image of 
the 3.5 cm continuum emission toward the center of the western clump, where a cometary 
shaped and several compact Hii regions, surrounded by a 0.1-pc halo, are seen. 

Since the minimum projected baseline lengths (0.78, 7.7, and 2.8 kA for the 3.5, 2 cm, 
and 3 mm data, respectively) and beam sizes (§2) are not identical, we reconstructed all 
the interferometric images with a I'.'S beam using visibilities with a uv-Yadius greater than 
7.7 kA. This makes our images insensitive to structures more extended than 27" (i.e. 0.46 
pc). Thus only the 0.1-pc halo common to all of the cluster members is resolved out. Figure 
2b represents an overlay of such maps at 3.5, 2 cm and 3 mm. Most of the 3 mm peaks 
are associated with previously identified cm sources (GRM85; WC89; Garay et al. 1998), 
suggesting that our 3 mm map is sensitive to both free-free emission from ionized gas inside 
the compact Hii regions and thermal emission from dust grains in the surrounding cores. 
To establish the properties of the dust continuum emission, we need to remove the free-free 
emission from the 3 mm map. We produced a "pure" free-free emission image at 3 mm 
by extrapolating the 3.5 cm map assuming optically thin emission [I^y oc z/^°'^). Then we 
subtracted this from the 3 mm continuum image, thus obtaining a "pure" dust emission 
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map (see Figure 2c). Taking into account the calibration errors at the various wavelengths, 
we estimate a subtraction uncertainty of 4.3 mJy beam~^. The resulting map presents a 
strong (29a) peak lying between the UC Hii regions labeled A and C (see the yellow cross in 
Figure 2a), where no emission peak is seen in 2 and 3.5 cm images. Instead, the 3 mm peak 
coincides with a molecular core (Figure 2d) detected in the CH3CH2CN (11-10) and SO2 
(83,5 — 92,8) lines: these molecular species are typical tracers of HCs (van Dishoeck & Blake 
1998). The core is also traced by other molecules such as NH3 (Garay et al. 1998), CH3CN 
(Kurtz et al. 2000), HCOOCH3, H^^CN, HN^^C, Hi3C0+, and SO (R. Furuya et al. in 
preparation). Such a richness in molecular species makes this core very similar to Orion-KL 
(Wright, Plambeck, & Wilner 1996) and other HCs (Kurtz et al. 2000). In addition, the 
presence of H2O (Hofner & Churchwcll 1996) and OH (GRM85) masers strongly suggest 
that the (proto)star deeply embedded in the core is too young to ionize its surrounding 
material because masers are believed to disappear during the development of an UC Hii 
region. No other prominent "dust" condensation was detected above a 5(7 upper limit of 
21.5 mJy beam^^, although extended emission is seen towards some of the Hii regions. Such 
an upper limit corresponds to a sensitivity in gas plus dust mass of 27 Mq beam"-^ for a 
typical HC temperature of 100 K (Kurtz et al. 2000). The large majority (~85%) of the 
HCs known to date (e.g. Table 1 of Kurtz et al. 2000) are above this hmit, but we cannot 
rule out the possibility to have missed the least massive ones (~10 Mq). 

We have identified the UC Hii regions in the cluster with a procedure, which iteratively 
fitted and subtracted prominent emission peaks from the cm maps, until the signal in the 
residual image was below 5a. The borders of the cometary shaped UC Hii region A and the 
source J at 2 cm were identified as the 5a contour levels enclosing them, whereas for all the 
other UC Hii regions a 2-D gaussian fit was used. This gave us the peak position, S^, and 
deconvolved source diameter of the sources. In the estimate of the error on S^,, we considered 
the formal error of the fit, the fiux calibration uncertainty, and the contribution due to the 
noise in the map. In this way, we identified 9 UC Hn regions, labeled A-D, F, G, I, J, and 
K. All of these were identified both at 3.5 and 2 cm, with the sole exception of J, too faint 
to be detected in the 2 cm image. 

In order to derive the properties of the UC Hii regions at 3 mm, it was necessary to 
disentangle the HC from UC Hii regions A and C. For this purpose, we used the free-free 
subtracted map at 3 mm (Figure 2c) to fit the HC with a 2-D Gaussian. Then this gaussian 
model was subtracted from the original 3 mm image. The resulting "clean" image was 
treated in the same way as the 3.5 and 2 cm maps to identify the counterparts of the cm 
sources. Because of the difficulty to define unique area for object F at 3 mm, we estimate 
an uncertainty equal to the difference between the flux estimate obtained from Gaussian 
fitting and that from integrating the emission inside the 5a contour. In Figures 2b and c. 
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one can see that an isolated 3 mm source lies to the east of A: we assume that this emission 

is associated with the tail of the cometary UC Hii region A rather than to a distinct YSO. 
We thus identified 6 objects at 3 mm, including the HC, above a detection threshold of 5o". 
In conclusion, we have identified a total of 10 high-mass YSOs: 1 HC and 9 UC Hii regions 
(Table 1). Note that we detected all the radio sources previously identified by other authors, 
except source E of WC89. All 9 UC Hii regions have been detected at more than 2 bands at 
cm wavelengths, if one takes into account identifications reported in the literature. 



4. Discussion 

We analyzed the continuum spectra of the whole region and of each cluster member 
by modeling the emission as originating from a dust-free Hii region surrounded by a dusty 
molecular shell. Spherical symmetry and homogeneity are assumed for both the Hii region 
and the shell. The spectrum of the whole region is shown in Figure 106 of WC89. To the 
data reported in this figure we have added our fiuxes at 3 and 1.2 mm (§3), and that at 
350 iim by Mueller et al. (2002). Most of the 3 mm emission is due to free- free, while the 
emission shortward of 1.2 mm is originating from dust. Given the complexity of the region 
seen in free-free emission, we estimate only a pc-scale mean dust temperature (Td) and a 
total mass of the gas plus dust (M^). Therefore, the contribution of the Hii region was 
considered only to subtract the free-free component from the spectrum. We adopted a dust 
absorption coefficient K,,y = ^^{I'/i'o)'^ with fi;o=0.005 cm^ g~^ at z/o=230 GHz (Preibish et 
al. 1993) and /3=1.5. We obtain rd~42 K, and Md~2800±100 M©, which are very similar 
to the values of other clumps containing high-mass YSOs (e.g. Fontani et al. 2002, and 
references therein). 

Subsequently we have analyzed the continuum spectra of the 10 sources in the cluster 
using the above model. In the fits we fixed radii of Hii regions (-Rhii) and dusty shells 
(i?d) to the values derived respectively from the 3.5 cm and 3 mm maps, and the electron 
temperatures (Te) to those obtained from recombination line observations (Garay et al. 
1998). For source A at the 3 bands, D and F at 3 mm, and J at 2 cm, we defined an effective 
radius as that of a circle with the same area as the UC Hii region or the dusty shell. We 
assumed a mean of 7200 K for the objects where Tc are not available in Table 5 of Garay 
et al. 1998. Therefore, free parameters were the emission measure {EM) of the Hii region, 
and Md and of the dusty molecular shell. Our spectral analysis shows that fluxes at 
3 mm (Table 1), 20.8, 18.1, and 11.7 /xm (De Buizer et al. 2003) towards objects B, F, I, 
and the HC are fitted with combinations of thermal dust plus free-free emission, implying 
that they are likely to be embedded in dust cocoons. Only the HC is not detected at cm 
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wavelengths, which indicates that the putative embedded stars (or star) are too young or not 
massive enough to develop an UC Hii region. We favor the idea that the HC contains very 
young massive (proto)star(s) as it is in all aspects (richness of the molecular composition, 
large mass of 800 M©, high temperature of 65 K) similar to the other HCs known in the 
literature, which are believed to be the cradle of early-type stars (Kurtz et al. 2000). If the 
HC contains a deeply embedded 0-B star, one may estimate the maximum radius of the 
associated Hii region from the sensitivity of the 3.5 cm image. Assuming that the free-free 
emission is optically thick, such a radius must be smaller than 2.4x10"^ pc. The continuum 
spectrum of D is rising from cm wavelengths to 3 mm (Table 1), but does not show MIR 
emission. This sets upper hmits on 7^ and M^. The 3 mm flux of source A (589±115 mJy) 
is only marginally above the value expected (409 mJy) from extrapolation of the free-free 
cm fluxes. Considering the uncertainties and the fact that its effective is smaller than 
that of i?Hii) we believe that thermal dust emission is not dominant at 3 mm. Finally, the 
spectrum of sources C, G, J, and K can be fitted with free-free emission only, suggesting 
that they have already dispersed most of the surrounding material. 

An unbiased survey such as ours could be used to estimate the lifetime of HCs with 
respect to that of UC Hii regions on the basis of number ratio between the two. However, 
this applies only to a series of star formation episodes randomly distributed in time, and 
not if one is observing the result of a coherent burst of star formation (e.g. Welch et al. 
1987). While in the former case the different numbers of sources in different evolutionary 
stages would be directly related to the time spent in each phase, in the latter hypothesis such 
numbers would be determined by the local physical conditions around individual stars. It 
is hence important to find out whether the ages of our sources are spread over a sufficiently 
large time interval or all very similar. 

Naively, one might expect that the age of an Hii region is strictly related to its diameter. 
However the correspondence between size and age is still matter of debate. Nevertheless, 
in most expansion models of an Hii region is related only to the spectral type of the star 
and to the properties of the surrounding environment. It is hence reasonable to assume that 
the expansion will come to an end when pressure equilibrium is reached between the ionized 
and the surrounding molecular gas. As illustrated by De Free, Goss, & Gaume (1998), their 
eqs. (4) and (7) give the final Hn region radius (-Rnif ) ^ ^ function of the stellar Lyman 
continuum photons, , and the density and temperature of the molecular gas. Here we 
derived spectral types of BO for A, B0.5 for B-D, F, G, I, K, and Bl for J. The value of R^f 
should be compared to the current UC Hn radius measured from the maps. Although it is 
impossible to obtain an "absolute" estimate of the age, ratios of Rmi/ Rmt be used 
to estimate "relative" ages of UC Hii regions, with lower ratios corresponding to younger 
objects. To compute R^^if, we used the corresponding values of and T^, and a mean 
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density derived from the SIMBA observations (~ 3 x 10^ cm~^). The ratio ranges from 0.3 
to 0.9 (see Table 1). If the UC Hii regions have not yet reached pressure equihbrium with 
the surrounding cloud, such a ratio seems to indicate that they are in similar evolutionary 
stages, thus supporting the "starburst" scenario. 

However, several caveats are in order: the density of molecular gas may vary over the 
natal cloud and the Lyman continuum of the stars may have been underestimated because of 
optically thick free-free emission and dust absorption inside the Hii region. This may easily 
introduce an uncertainty of an order of magnitude on R^lf thus making the "starburst" 
hypothesis questionable. In fact, the expansion timescales obtained from cq. (3) of Dc Prcc, 
Goss, & Gaume (1998) span a small range (120-480 yr), suggesting that the UC Hn regions 
are close to pressure equilibrium. It is hence worth considering also the opposite scenario, 
and take the number ratio between HCs and UC Hii regions as an estimate of their relative 
lifetimes. 

Considering Poisson statistical errors, the number of HCs and UC Hii regions are re- 
spectively 1±1 and 9±3, so that their ratio is This imphes that the HC phase 
(thc) should last less than ~l/3 of the UC Hii region phase (thu), with a probabihty of 4% 
that the ratio exceeds this value. This conclusion seems to differ from that inferred for the 
0-star cluster in W49N (Wilner et al. 2001), where rHc/rnii - 25-100%. The two results 
could indicate that thu — 3-4 thc- Alternatively, the difference might be easily explained 
in the "starburst" scenario, where the ratio between HCs and UC Hii regions depends on 
thc environment around each star. Finally, another caveat is in order. Our statistics might 
be affected by the limited sensitivity of our 3 mm map. As previously explained, we are 
sensitive to HCs with masses in excess of ~27 Mq, whereas at cm wavelengths we detect 
stars earlier than Bl. If O stars are born in massive HCs and B stars in lower mass HCs, 
our comparison between UC Hii regions around B stars and HCs hiding O stars is bound to 
be unreliable. More numerous and sensitive interferometric observations of 0-B star clusters 
will be needed to improve the statistics. 

The authors gratefully acknowledge an anonymous referee whose comments significantly 
improved the quality of the paper. R. S. F. thanks Prof. A. I. Sargent and J. M. Carpenter 
for critical readings of the manuscript and encouragement, M. Niclbock and L. Haikala for 
their generous help at SIMBA observations and data reduction. Many thanks arc due also 
to T. Saito for his early contribution to this study, and to thc staff at OVRO, NRO, NRAO, 
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Science Foundation through NSF grant AST 02-28955. 
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Fig. 1. — Parsec scale view of the cloud hosting the G19.61-0.23 region in 1.2 mm continuum 
emission. Contours start from the 3a level of noise level (30.7 mJy beam~^ ), and end at 
the 15(7 level with 3a steps. The hatched ellipse in the bottom right corner indicates the 
SIMBA beam size (22"), and the central white box marks the region shown in Figure 2. 

Fig. 2. — (a) 3.5 cm continuum emission map taken with VLA. The hatched ellipse in 
the bottom right corner indicates the beam size. The white contour corresponds to the 
10a level. The yellow cross indicates the position of the dust continuum emission peak 
(R.A.= 18^27'^38.'066, Decl.= -11° 56' 37'/20 in J2000). (b) Overlay of 3 mm continuum 
emission map taken with OVRO and NMA (grey scale plus black contour) on the 3.5 cm 
(thick purple contours) and 2 cm (thin green contours) maps CLEANed with the 3 mm 
beam (see text). The labels denote identified sources. Contour levels for the 3 mm map are 
5a intervals starting from the 5a level ((Jamm = 4.3 mJy beam~^ ). For the sake of clarity, we 
plot only 5a, 10a, and 30cr levels for the 3.5 cm map, and 5a, 10a, 20a, and 40a levels for 
the 2 cm map, where a =0.38 and 1.6 mJy beam^^ , respectively, (c) "Pure" dust continuum 
emission map obtained after subtraction from the 3 mm image of the 3.5 cm map suitably 
scaled to 3 mm (see text). Contours start from the 5a level with a 5a step {adust = 4.3 mJy 
beam~^ ). (d) Overlay of total integrated intensity maps of CH3CH2CN (11-10) (grey scale) 
and SO2 83^5 — 92,8 emission (green contours) taken with the NMA. Both the grey scale and 
green contours start from the 5a level with 5a steps, where noise levels and beam sizes are 
14.1 mJy beam-^ and 6'.'40x3'.'65 for CH3CH2CN, and 9.4 mJy beam'^ and 4'.'59x2'.'95 for 
SO2. The open red circles and 6 blue triangles indicate the positions of the OH (GRM85), 
and H2O (Hofner and Churchwell 1996) maser spots, respectively. The yellow cross has the 
same meaning as in (a). 
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Table 1. Parameters of Identified Sources 
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^From the left, source labels identified in this work, WC89, Garay et al.(1998), and De Buizer et al. 
(2003). 

''Flux densities are measured with a 1'.'5 beam and a minimum projected baseline length of 7.7 kA (§3), 
and those in parenthesis are marginal detections that the significance levels of the peak intensities range 
between > 3cr and < 5cr. All the results in the following analysis based on the marginal detections are 
shown in parenthesis. Upper limits are calculated with the lOcr level noise levels over the 1'.'5 beam. 

■^Measured at 14.9 GHz for the sources I, G, and K. 

dEffective radius (§4) 

''Upper limit is from lack of well-defined emission peak at the 18.1 /um map (Figure 8 of De Buizer et 
al. 2003) 

^Assuming optically thick free-free emission (§4) 

^Obtained from the 18.1 /im map (Figure 8 of De Buizer et al. 2003) 

''The source J is located inside a circle given by the position of the source CD in G98 and its angular 
size (see their Table 2). The source K corresponds to the north-western peak of the 6 cm source D of 
WC89. 



